Twirling DNA Rings - Swimming Nanomotors Ready for a Kickstart 
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We propose a rotary DNA nanomachine that shows a continuous rotation with a frequency of 
10 2 — 10 4 Hz. This motor consists of a DNA ring whose elastic features are tuned such that it can 
be externally driven via a periodic temperature change. As a result the ring propels itself through 
the fluid with a speed up to microns per second. 

PACS numbers: 05.40.-a, 47.15. Gf, 87.15. He 



The long lasting dream of scaling mechanical devices 
and machines down to the nanoscale (as popularized 
by R. Feynman [if and carried on by several visionary 
groups worldwide |2|]) continues to fire the imagination 
of researchers - now in the third generation. Among 
many experimental difficulties that appear in this con- 
text, choosing the proper material for the assembly of a 
nanodevice turns out to be crucial. Important material 
requirements are: stability, self-assembly ability, modu- 
larity, replicability, switchability, experimental tractabil- 
ity. Presently one of the most promising materials ful- 
filling those requirements is DNA |3|. Assemblies based 
on DNA hybridisation chemistry_0, 0, S as weu as 
conformational DNA transitions [8j were successfully ex- 
ploited to generate periodically switchable nanodevices. 
However, despite their beauty and conceptual original- 
ity all of these devices suffer one major problem: the 
large kinetic barriers involved in the switching process 
boost their switching time per cycle to ~ 10 3 s, four or- 
ders of magnitude slower than their natural counterparts 
(biological molecular motors). A natural question arises 
then: Can one achieve subsecond switching times with a 
DNA nanodevice? Can such a device be operated in some 
manner to swim directionally and faster than 1/im/s? In 
this letter we show theoretically the principal feasibility 
of such DNA nanomachines. 

Let us in the following propose a surprisingly simple 
nanomotor: a DNA miniplasmid, cf. Fig. 1(a). We will 
show that despite its structural simplicity a miniplasmid 
can be turned into a nanomachine able to produce fN 
forces and self-propelling at speeds of several microns per 
second through the fluid. In order to run the plasmid as 
a motor we use here the Euler-angle %jj (cf. Fig. 1(a)) as 
the relevant degree of freedom [9j. The main idea now 
is to induce a directed current (ip) - in a manner sim- 
ilar to the rotation of a closed rubber tube around its 
central circular axis - via non-equilibrium fluctuations 
and the ratchet effect [3 EH, cf. Fig. 1(b). As a re- 
sult the twirling ring generates a hydrodynamic flowfield 
(shown in Fig. 1(c)) that remarkably induces an efficient 
self-propulsion of the motor as detailed below. 

The elastic distortion energy of a DNA ring with ra- 
dius R parametrized by the arc length parameter s will in 
general be described by three Euler angles 9 (s) ,<j> (s) and 




FIG. 1: The operation principle of the DNA-minicircle pro- 
peller: (a) The twirling degree of freedom, (b) The elastic 
energy as a ratchet potential, (c) The flowfield around the 
twirling ring induces its translational velocity v z . 



ip (s) via E e i = \k B T f Q 71 £)i=i,2,3 h - Hif ds with 
u>i = 4>' sin sin t/j + 8' cos 4>, uj?, = cf>' sin 8 cos if) — 0' sin ip 
and 0J3 = (j)' cos# + ip' [12|]. Here l\ and I2 are the 
two principal bending persistence lengths and k\ and K2 
are intrinsic curvatures in two corresponding perpendic- 
ular directions. I3 denotes the twist persistence length 
and - for simplicity - we choose K3 = 0. The parame- 
ters Ki and U reflect the anisotropic bendability as well 
as intrinsic bendedness of the plasmid; here for sim- 
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FIG. 2: The rotational current {if>) and the induced transla- 
tional velocity v z as a function of the dimensionless frequency 
/ of the temperature (potential) oscillations. The DNA ring 
has the following parameters: R — lOnm, ro = lnm (typi- 
cal DNA minicircle), Zi = 45 nm, I2 = 50 nm, Ki = k,2 — 
(200nm) _1 leading to the ratchet potential displayed in the in- 
set. Displayed are the asymptotic expressions, Eqs.|5] (dashed 
line) and|S| (dashed-dotted line) together with the numerical 
solution of Eq. [5] (thin line) for a temperature ratchet with 
At = 0.03. The thick solid line corresponds to an oscillating 
potential ratchet with Ae = 0.3. See text for details. 

plicity we assume them to be independent of the arc- 
length throughout the molecule (e.g. for a 10 basepair 
(bp) periodic sequence). For the case of DNA mini- 
circles of short length (2ttR < k) and with constant 
Ki and li fulfilling the weak bending anisotropy condi- 
tion max{|Zi — Z2I /R, iiKij I2K2} *C I3/R only the con- 
formations close to the circular untwisted state will con- 
tribute, i.e., those close to 8 (s) = tt/2, 4>(s) = s/R and 
ip (s) = const. This leads then to the required ratchet 
potential acting on tp: 

— 7-^7- = — cos (2ip) + 2liKi costp — 2I2K2 sin?/; (1) 

ttkbT 2R 

From Eq.^we see that for generating a ratchet potential 
we need both nonzero bending anisotropy, l\ — I2 7^ 0, 
as well as non-vanishing intrinsic curvatures, Ki. 2 7^ 0. 
The inset in Fig. 2 demonstrates that reasonable small 
values of anisotropy and intrinsic curvature can induce a 
well-defined ratchet potential. 

The Fokker-Planck equation describing the time evolu- 
tion of the probability density P (ip, t) of the Euler angle 
ip writes 

BP d fdE el , dP\ 

with the friction constant ( that we will compute below. 
As a source of non-equilibrium we will choose here a time- 
dependent variation of temperature T (t) , cf. Ref. 0] . 

Before we compute the friction constant C we need to 
shed some light on the low Reynolds number hydrody- 
namics of the twirling DNA ring. The latter turns out to 



be peculiarly related to the inviscid (ideal) fluid vortices 
(rings of smoke) and as a matter of fact both of them 
propagate in analogous manner. To see this we first re- 
mark that for a reasonable ring radius R = 10 nm (a 
typical miniplasmid of w 200 bp) and the DNA helix ra- 
dius ro = 1 nm the slender body approximation [l4| is 
valid with the slcnderness paramentcr e = tq/R = 0.1. 
In the spirit of the slender body theory one approximates 
the flow-field around the twirling ring by superimposing 
rotlets u rot (x; s) = x (x - c (a)) / |c (s) - x| 3 

placed along the ring centerline c (s) with arclength pa- 
rameter s. The rotlet strength T = ^w c r o ^ s given 
in terms of the angular velocity lo c of the ring about 
c(s). The full velocity profile is then given by u (x) = 
J^ R u rot (x;s)ds; cf. also the stream lines around the 
rotating ring shown in Fig. 1(c). When integrating u (x) 
over the DNA ring (slender torus) surface in the limit of 
small ro / R one obtains a net translational velocity in the 
z-direction: 

The fact that Eq. |31 coincides with the well known ex- 
pression from ideal flow vortex theory 0] should not 
surprise if we recall that a rotlet u rot (x; s) is nothing 
else but the expression for the velocity field of an ideal 
point vortex. But despite this kinematic analogy between 
the twirling DNA and an ideal vortex ring, dynamically 
they are quite different. The propagation of an ideal 
vortex ring does not require any external forces/torques 
and is governed by conservation of kinetic energy and 
momentum. In sharp contrast to that the low Reynolds- 
number (Stokes) flow is governed by dissipation and the 
motion of twirling DNA ring requires the action of a 
torque A^ c = 8n 2 XQriRuj c (77 = 10~ 3 Pa s, the water 
viscosity) about the central axis c. The latter expres- 
sion can be verified by integrating the tangential stresses 
generated by u (x) over the ring surface. More generally 
by virtue of the linearity of the Stokes equations we can 
derive a resistance matrix (My) relating the angular ve- 
locity lj c (about the circular axis c) and velocity v z (in 
the z-direction) with the corresponding external torque 
N c and force F z : 

( F z \ 2 / Mu Mia \(v s \ 

Combining the previous expressions obtained for v z (w c ) 
and A^ c (ui c ) (F z = 0) together with the result of 
Johnson and Wu [l^ for the drag on a rigid slender 
torus we obtain entries in the leading order: Mn = 
2i?(ln8/e + 1/2) -1 , M 22 = 2r%R and M 12 = M 2 \ = 
rl (ln8/e - 1/2) (ln8/e + l/2)~\ Note that the symme- 
try of the resistance matrix being a general feature of 
swimmers in the Stokes flow |l7j provides a good check 
for the consistency of the involved calculations. Having 
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the mobility relation 0] one can consider different types 
of motion, e.g.: (i) The ring is twirling freely (N c = 0) 
and moving under the force F z (or at fixed v z ). (ii) The 
ring is prevented from rotation (u) c ~ 0) and moved by 
the force F z (or at fixed v z ). (Hi) The ring is held in po- 
sition (v z = 0) by imposing a force F z counterbalancing 
the action of torque N c . (iv) The DNA ring is free to 
move (no external force applied, F z = 0) under an exter- 
nally imposed torque N c (or at given lo c ). Comparing, 
for instance, the velocities v z and v z from cases (i) 
and (ii) one sees that v z /«* — 1 = Mj^/detM/y ~ 
(e/2) 2 In (8/e); i.e., a ring with an isotropic DNA se- 
quence (able of twirling) settles faster than a ring with 
very high barriers to twirling (uj c = 0), though in practice 
the difference is negligible, for instance ~ 1% for e = 0.1. 
By comparing cases (Hi) and (iv) we conclude that a ring 
twirling at fixed lu c and forced not to translate requires 
a slightly larger torque N™ 1 (w c ) than the unconstrained 
freely translating ring with N™ (ui c ). However the rela- 
tive difference is again small and of order O (e 2 In (8/e)J . 
Therefore by dropping this marginal correction to the 
leading order we obtain in both cases (Hi) and (iv) the 
angular friction constant ( = N c (oj c ) /u> c » S^-qr^R, 
the quantity that appeared above in Eq. [21 Note that 
the latter is the same (in our e <C 1 leading order ex- 
pansion) as for a straight cylinder with radius tq and 
length 2irR. Finally another interesting feature that can 
be read off Eq. 0] is the efficiency of the twirling ring 
propulsion. The latter is independent of the mechanism 
of twirling and can be defined as the ratio of the power 
Po = 2ir 2 r]Mi 1 v1 dissipated by a (for simplicity) rigid 
ring directly moved by a force as compared to the power 
P twir i = ^N c ui c dissipated by twirling propulsion at the 
same translational speed. For a ring with R = 10 nm 
we have Po/Ptwirl ~ 0.8% which is comparable to the 
efficiency of bacterial propulsion by a rotating flagellum 

Having determined the friction constant £ we return to 
the ratchet dynamics given by Eq. [21 with the twirling po- 
tential Eq. ^ To obtain the directed twirling freque ncy 
lu c := W = -§(9§pP + k B T^) we follow Ref. Q 
by choosing a periodic time dependent temperature vari- 
ation as follows: T (t) = T [1 + A T aia(2n frt)] with T 
the mean temperature, At the relative amplitude and 
fx the frequency of the temperature oscillation. For 
the case of fr sufficiently larger than the inverse of the 
characteristic relaxation time tq = 47r 2 £/ (ksTo) of the 
twirling degree of freedom (but still much smaller than 
the frequency of average thermal molecular kicks) an 
1//t asymptotic expansion for the current (ip) can be 
employed (ljj . After a long calculation we obtain (for 
/t > /res) (ip) up to terms of order O (/-p 3 ): 

/ . \ 12n 3 A 2 T (k B T f (l 2 - h) [R 

\W EeiW n„ Beit*) ^ > 



From Eq. [S] we see that for an isotropically bendable 
DNA sequence (Z2 = ^1) the directed current vanishes. 
The same is true if the intrinsic curvature direction co- 
incides with one of the principal axes (i.e. if K\ or K% 
vanish). Both observations are intuitive as in either case 
the ratchet potential, Eq.^ becomes left-right symmetric 
and the ratchet effect disappears. 

The low frequency adiabatic limit is obtained from the 
asymptotic expansion of P(tp,t) (-0 = *P/2n, i = t/ro) 
for small parameter / = /tto, i.e., P«Fo + fPi + f 2 P2- 
Rather involved calculations lead to 

U\ =Jlf dt^FWi (6) 
\ / t Jo F 

with E$,t) = F$,t)- 1 = e -E c iW/k B T(t) and the ab _ 

breviations (...) and (...) defined as in 0] but with the 
integrations with respect to "0. Furthermore the den- 
sity distributions Po and Pi from the upper expansion 
are given by Pq = E / E (Boltzmann distribution in the 

adiabatic limit) and Pi = ^-E (FcI - iPPct) with 

Equations [S] and together with Eq. ^ an d C = 
87r 2 r]rQR allows us to get the twirling speed uo c = (ip) 
and by virtue of Eq. [3] the induced translational velocity 
v z (lu c ) for arbitrary DNA elastic parameters ^=1,2 and 

^i=l,2- 

How fast can we operate the twirling ring machine? 
We shall assume some realistic parameter values for the 
DNA ring: R = 10 nm, tq = 1 nm (typical DNA mini- 
circle) leading to ( = 2 ■ 10~ 7 fcsTs. Furthermore we set 
li = 45 nm, I2 = 50 nm, K\ = K2 = (200nm) _1 which 
corresponds to a rather modest anisotropy and intrinsic 
curvature. For the temperature variation amplitude we 
choose AT = ±10 K, i.e., At ~ 1/30 (at room temper- 
ature To = 300 K). Figure 2 provides a log- log plot of 
the rotational current and the corresponding drift speed 
of the ring as a function of the dimensionless frequency 
/ of the temperature variation. The thin solid curve 
gives the numerical result obtained from Eq. [21 the two 
straight lines correspond to the analytical results for the 
two asymptotic cases, Eqs.EJandEl As can be seen from 
this plot the two limits show a f~ 2 and f 2 dependence, 
respectively, in accordance with Eqs.[31and[S| The maxi- 
mal rotational current is achieved in the crossover region, 
namely uj c ~ 200rad/s for / sa 10 _1 . Following Eq.^lthis 
implies a translational velocity of v z = 50nm/s. 

Such fast temperature oscillations are technically fea- 
sible and might be most conveniently generated by adi- 
abatic pressure variations e.g. by ultrasound. An- 
other promising method is to use the inductive heat- 
ing of metal nanocrystals that are covalently attached 
to the DNA ring. In fact, this method has been suc- 
cessfully used to control the hybridization behavior of 
DNA |20(. This might also point towards an alterna- 
tive way of driving the ratchet, namely via a periodic 
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variation of the elastic properties of the ring. Operating 
the system close to the DNA duplex melting tempera- 
ture is likely to induce strong oscillations in the overall 
ring stiffness. The thick solid line in Fig. 2 shows the 
rotational current obtained when the elastic energy is 
varied as E e i (tp, t) — E e i (tp) (1 + Ae sin (27t/e£)) where 
we chose the relative amplitude Ae = 0.3. As can be 
seen from Fig. 2 the maximal current of this oscillat- 
ing potential ratchet occurs roughly at the same fre- 
quency as that of the thermal ratchet but the value of 
u> c is much higher, namely on the order of 2 x 10 4 rad/s 
which implies a quite notable translational velocity of 
v z = 5/im/s. As a comparison a typical bacterium moves 
at 30/im/s. Our ring ratchet (with oscillating potential) 
resembles in many respects "real" biological nanomotors. 
Besides its nanoscopic size (radius 10 nm), swimming ef- 
ficiency (0.8%) and speed (4^im/s) it can generate forces 
and torques close to that of biomolecular motors. Al- 
though the net translational force resulting from Eq. 0] 
F z = 4tt 2 t]Mi2uj c w 0.6 fN is comparably small (due 
to cancelling of most of the stresses), the local torque 
iV c = 87r 2 ryrQ_Rti> c w 0.004/cbT and the force acting at 
the DNA surface F\ oc = N c /r — C,uo c /r a » 16 fN are 
significant if we consider the simplicity of the mechanism 
behind. 

From an experimental point of view one should be 
aware of the fact that a ring (twirling or non-twirling) 
looses its initial orientation almost instantaneously due 
to rotational diffusion. The typical relaxation time scale 
of this process is on the order ijR 3 / (fceT) (up to loga- 
rithmic corrections 21]) which for a ring with R = 10 nm 
leads to 10~ 7 s. That means that a single twirling ring 
in solution will not perform any noticeable translational 
drift. A possible solution to the problem is to put the ring 
on a "track", e.g. to thread it on a straight DNA chain. 
Another promising direction is to prepare a semi-dilute 
solution of such rings and then study their response due 
to an induced twirling. The resulting hydrodynamic flow 
can lead to an attraction between the rings that might 
facilitate interesting collective behavior of the nanomo- 
tors like the formation of columns of twirling rings that 
drive a solvent flow through the resulting channel. 
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